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Abstract 


A  version  of  statistical  theory  is  applied  to  the  hydrogen 
molecule-ion,  with  good  results  for  the  electronic  energies  of 
the  lowest  attractive  and  repulsive  states.  For  the  former,  a 
binding  energy  of  0.085  a.u.  is  obtained  with  an  internuclear 
distance  of  2.10  a.u.,  to  be  compared  with  0.103  a.u.  and  2.00 
a.u.  obtained  from  the  exact  solution  of  Schroedinger's  equation. 
The  method  employed  involves  the  use  of  hydrogen  atomic  orbitals, 
although  when  the  latter  refer  to  unbound  states  they  are  approxi¬ 
mated  by  plane  waves. 
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1.  Introduction 

(l  2) 

The  first  application  of  the  Thomas-Fermi  theory  *  to  molecular  systems 

(3) 

occurred  soon  after  its  introduction,  when  Hund'  used  it  to  approximate  the 
electron-density  and  the  electric  field  in  the  nitrogen  and  fluorine  molecules. 
Since  then,  the  original  statistical  theory  of  Thomas  and  Fermi  has  been  em¬ 
ployed  to  estimate  the  binding  energy  of  molecular  systems^.  The  numerical 
results  which  have  been  obtained  from  the  use  of  the  original  Thomas-Fermi 
theory  have  been,  at  best,  only  in  fair  agreement  with  experimental  values  of 
bond  energies  and  bond  distances;  at  worst,  an  absence  of  binding  has  been  in¬ 
dicated.  Recently,  the  limitations  of  the  Thomas-Fermi  theory  in  this  connec¬ 
ts 

tion  have  been  established  by  Teller  .  He  has  shown  that  this  theory  is 

incapable,  in  terms  of  its  original  formulation,  of  yielding  energetically 

stable  molecular  systems.  A  similar  conclusion  is  obtained  for  the  version 

of  statistical  theory  known  as  Thomas-Fermi -Dirac  theory,  in  which  the  effect 

of  electron  exchange  has  been  incorporated^  into  the  original  theory. 

In  the  years  following  its  introduction,  the  theory  of  Thomas  and  Fermi 

(7) 

has  been  subjected  to  a  variety  of  other  modifications  .  Of  all  of  these, 

/  D\ 

the  one  due  to  von  Weizsacker'  '  has  produced  what  seems  to  provide  the  most 
radical  and  most  effective  connection  of  the  original  theory  with  quantum 
mechanics^.  With  this  modification,  Goinbls^^  has  obtained  approximate 
values  for  the  binding  energy  and  the  bond  length  of  the  nitrogen  molecule 
which  are  in  good  agreement  with  the  experimental  values  of  these  quantities. 
From  this  result,  a  sufficiently  quantum-mechanical  version  of  statistical 
theory  appears  to  be  essential  for  a  realistic  application  to  molecular  sys¬ 
tems. 

An  alternative  to  the  von  Weizsacker  theory  is  one  which  is  the  formal 
analogue  of  the  Thomas-Fermi  theory,  but  employs  a  basis  of  eigenfunctions  of 
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some  suitable  approximate  Hamiltonian  which  is  more  closely  related  to  the 
physical  system  of  interest  rather  than  the  basis  of  eigenfunctions  of  momentum 

f 

implicit  in  the  original  theory.  Such  a  version  has  been  shown  to  yield  rela¬ 
tively  good  results  for  the  energy  values  and  electron  densities  of  helium-like 
atoms The  present  paper  is  concerned  with  an  application  of  this  version 
of  statistical  theory  to  approximate  the  binding  energy  of  the  simplest  mole¬ 
cular  system  —  the  hydrogen  molecule-ion.  Because  of  the  simplicity  of  this 
system,  the  usual  complications  attending  the  determination  of  a  self-consistent 
field  can  be  avoided.  As  a  consequence,  the  inherent  limitations  of  the  theory 
can  be  best  assessed  under  the  most  stringent  circumstances. 


2 .  Formalism 


In  the  interest  of  brevity,  we  shall  simply  sketch  the  formalism  to  be 

(l2) 

employed  herev  .  The  properties  of  a  single-particle  (i.e.,  an  electron, 
since  we  tacitly  suppose  nuclei  of  infinite  mass)  can  be  determined  from  the 
density  matrices  (representatives)  of  the  system  determined  by 

oO 

?M(x',  x)  =  21  ”tn*^’  ^  ^  ^  M  "  5)  »  (2,l) 

n*l 


where  x  stands  for  the  configurational  and  spin  coordinates  of  the  particle, 

H  is  the  Hamiltonian  of  the  system,  \^>n(x)'\  is  any  complete,  orthonormal  set 
of  functions  with  the  correct  boundary  conditions.  Here,  T^(y)  is  the 
Heaviside  unit  function  defined  by 

^(y)  =1,  y>  0, 

(2.2) 

-0,  y  <  0; 

a  convenient  representation  of  the  ^-operator  —  the  spectral  operator  of  H  — 


is 


V+ieo 

US  1  p  C  dz  M  -zH  n. 

m  -  JP  =  m  P  J  T  6  e  Y  >  0-, 

/  -i°2> 


the  parameter  Xm>  which  plays  the  role  of  a  "chemical  potential  per  state 


is  chosen  so  that 


Tr(^j,j)  =  \  dx  j^Cx,  x)  =  M,  a  positive  integer.  (2.I4) 

(13) 

With  this  normalization,  M  refers  to  the  M  lowest  eigenstates  of  H. 

In  terms  of  the  eigenfunctions  of  H,  (*„<*>).  the  density  matrices  assume 
the  form  originally  introduced  by  Dirac ^ 


a  -  2  -£<>.„  -  En>  0n*<a)  »)„(*) 


0  *(x')  0  (x), 


so  that  for  any  arbitrary  observable  A,  the  quantity 


<*>„  ■  ^  ^  I A 1 0n > 

/v-  n*=l 


yields  the  sum  of  the  expectation  values  of  A  for  the  M  lowest  eigenstates  of 
H.  In  particular,  we  see  that 


J.J. 

<E>m  -  Tr(fMH)  -  2  En, 


n=l 


where  —  are  the  eigenvalues  of  the  Hamiltonian.  It  can  be 

verified  that  Eq.  (2.6)  is  obtainable  fron/1^ 

M  M 

<*>„  -  j  X(»  dM-  J  A{Trfl,jd{Tr|M} 


(2.7) 
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which  emphasizes  the  role  played  by  A  ,  noted  earlier. 

Equations  (2.1|)  and  (2«7)  involve  the  traces  of  the  density  matrices  and 
are  thus  invariant  to  the  choice  of  basis  employed  for  this  evaluation.  They 
are  exact.  Moreover,  in  the  approximation  dealt  with  here,  they  will  be  used 
to  provide  an  estimate  of  the  electronic  energy  of  the  system  of  interest. 

Before  considering  any  approximation,  however,  we  examine  the  effect  of 
the  symmetry  of  the  system  upon  its  density  matrices.  Since  the  Hamiltonian 
of  the  system  is  left  invariant  as  a  consequence  of  the  application  of  certain 
symmetry  operations,  the  density  matrices,  by  Eq.  (2.1),  likewise  must  remain 
invariant  to  their  application.  In  particular,  the  inversion  symmetry  of  the 
hydrogen  molecule-ion  will  be  exploited  to  simplify  the  density  matrices  for 
that  system.  In  this  case,  we  have  the  Hamiltonian,  in  atomic  units,  of 
electronic  motion  (see  Fig.  l) 


H  * 

A< 


(2.8) 


Clearly,  the  Hamiltonian  remains  invariant  as  a  result  of  the  application  of 
the  unitary  transformation  $  ---  where 


-r  , 
o  ' 


(2.9) 


corresponding  to  inversion  through  the  mid-point  on  the  internuclear  axis. 
Furthermore,  since 


(2.10) 


the  inversion  operator  also  is  Hermitian.  Hence  its  eigenfunctions  comprise 
two  distinct  symmetry  speciess  those  which  are  symmetric  (gerade)  and  those 
which  are  antisymmetric  (ungerade)  with  respect  to  the  inversion  operation. 
Because and  H  commute,  their  simultaneous  eigenfunctions  possess  the  same 

r-  /Ni 

symmetry . 
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In  terms  of  such  simultaneous  eigenfunctions,  we  may  construct  two  sorts 
of  density  matrices.  Their  construction  is  facilitated  by  the  introduction  of 
the  orthogonal  projections 

s(±)  =  i(i±^).  (2.11) 

Since  and  H  commute,  we  may  obtain 

and,  ultimately, 

fM(x',  x)  -  x)  +  fM(-)(x',  x),  (2.12) 

where,  by  Eqs.  (2.1)  and  (2.9)  and  the  introduction  of  the  coordinates  of 
Fig.  1, 

j)  -  j  2  „<-?„>}• 

^  n=i 

(2.13) 

The  two  classes  of  density  matrices  and  are  demonstrably  orthogonal 

and  refer  to  the  inversion-symmetric  and  inversion-antisymmetric  states,  re¬ 
spectively,  of  the  hydrogen  molecule-ion.  It  is  noteworthy  that  the  resulting 

(l^) 

basis  functions  here  are  of  the  Pauling  type  ,  but  do  not  comprise  an 
orthonormal  basis. 

Analogous  to  Eqs.  ( 2 .14)  and  (2.7),  we  obtain 
M(i)  -  Tr(  ?„<*>>  -  i  2Z  <+n(?0)|-9'(X„<£)  -  H)  lin(?„)/’ 

r-j  nml 


t  *  2  -  ®|  +„<?„)>• 


(2.1b) 
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<E>„1')  -  Ir(fH(i)  J.J2  <l|»nff0)  I  H-S-a/)- H)  |  4n(?0)> 


n**l 

oo 


*.  i  2  <!(-?„>  I  B-3(XM(t)-H)|  +„(?„)>•  (2.15) 


3.  Quasi-classical  Approximation  —  Discrete  States 


The  approximation  to  be  considered  here  can  be  expressed  simply  as 


e 


-zH 

A/ 


-zV  -zHq 
e  e  ~ 


(3.1) 


where 


1 


(3.2) 


and 


(3.3) 


This  warrants  taking  \y  n(rQ) J  to  be  the  complete  orthonormal  basis  of 
eigenfunctions  of  the  hydrogen  atom  corresponding  to  one  of  the  nuclei  of  the 
hydrogen  molecule-ion.  With  this  identification,  we  find  it  convenient  to 
refer  to  nucleus  -2  of  Fig.  1  as  the  origin  and  obtain 


(3  JO 


where  E^0  is  an  appropriate  eigenvalue  of  the  hydrogen  atom.  A  corresponding 
normalized  eigenfunction  is^^ 


‘  Vr2>  ®a.  (92> 


and 


(3.5) 


o 


(3.6) 
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(The  azimuthal  angle  0^  measures  rotation  of  the  electron  about  the  inter- 
nuclear  axis,  but  is  not  included  in  Fig.  1.)  The  effect  of  inversion  of 
the  eigenfunctions  can  be  seen  to  be 

£  t(?2>  ■'  Vr3>  <e3>  •i’(02*’'V75r 

■  <-1)”  W  ®£m(93> 


‘  fn(-?3>- 

Combining  these  results  with  Eq.  (2.3),  we  obtain 

-  $  W  ;  *  75  *  ?:> 

2n  j 


(3.7) 


(3.8) 


The  approximations  to  Eqs.  (2. lit)  and  (2.15)  can  now  be  effected,  'kith  Eq. 
(3.7),  we  obtain 

Tr(f„(!))  -  i  X  2  P  -h  *  jr>l+n<?2>> 

a*  n=l  Jt=o  w  2n  3 

*  75  *  7>l+n(?2)>}  •  «•’> 

2n  i 

(17) 

Making  use  of  the  addition  theorem  for  spherical  harmonics  ,  one  can  obtain 

+1 


and 


(£+  i)  P£  (-cos  a)  -  2^  ®Aln(e2)  ®iim(e3)  e  }  * 

m»-X 


With  these  relations  it  can  be  established  that  Eq.  ( 2 .lit )  is  approximated  by 
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“d^  ^  jS  {Qj*  i){^Rnll(r2)l^(^M(")  +  ^2  +  7^\  Rn£(r2)/> 

t  (-if  <^RnA(r3)|  P^Ccos  a)-^(XM(±)  +"^2  +7-)|  Rn4(r2)/>}  ,  (3.10) 


2n  3 

2 

where  use  has  been  made  of  the  relation  ^  (-x)  ■  (-1)  ^  (x).  In  equivalent 

terms,  Eq.  (2.13)  is  approximated  by 

n-1 

S 

Z  =o 


<■•>;•’  -  -  i  5  s  m y  v-,>> 

1  ("l)XRn4(r3)l  Vcos  a)^(  V_)  +  ^2  +  FJ}i  Vr2}>} 

-f,  |  S  <*♦  ={(WI  W1’  ^4)lw> 


*  «>  ^  ■ 

(3.11) 


Because  of  Eq.  (2.2),  the  integrands  in  Eqs.  (3.10)  and  (3.1l)  have  non-zero 
values  in  restricted  regions  of  configuration  space.  In  particular,  such  non¬ 
zero  values  occur  for 


Hence,  defining 


k«(t)  --  -*.c±>  >  °> 


(3.12) 


we  must  restrict  the  various  integrations  to  the  region 


0  4  r3  ^  L""W  J~  * 

"  2n2 


(3.13) 
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It.  Approximation  for  the  Continuum 

The  expressions  of  the  previous  Section  are  formally  incomplete  transcrip¬ 
tions  of  the  results  of  Section  2.  They  include  only  the  eigenfunctions  of 
negative  energy  of  the  hydrogen  atom.  To  complete  the  transcription  one  re¬ 
quires  the  inclusion  of  the  eigenfunctions  of  the  continuum  of  positive  energy 

/  O  \ 

values.  Although  these  are  well-known'  ,  they  are  extremely  difficult  to 

handle  in  the  present  context.  As  a  consequence,  an  additional  approximation 

will  be  made  here  in  which  the  plane-wave  basis  of  the  original  Thomas-Fermi 

theory  will  be  used  for  the  continuum.  Such  usage  can  be  justified  only  in  an 

asymptotic  (E°  - *■  +  °°)  sense,  however.  Even  so,  the  use  of  the  original 

Thomas-Fermi  density  to  approximate  the  continuum  is  bound  to  introduce  errors 

related  to  the  poor  values  of  the  binding  energy  it  yields  for  atoms.  For  this 

(a) 

reason,  we  have  adapted  a  version  of  statistical  theory  '  which  has  been  shown 

to  yield  improved  results  for  atoms.  In  this  theory,  the  density  matrix  of  the 

system  is  developed  in  a  power  series  in  Planck’s  constant  and  terms  up  to  and 
2 

including  h  are  retained. 

In  essence,  the  adaptation  consists  in  representing  the  unbound  states  of 
the  hydrogen  atom  centered  at  nucleus  -2  by  eigenfunctions  of  momentum  and 
their  corresponding  energy  by  their  (asymptotic)  kinetic  energy.  Compared  to 
Eqs.  (3.5)  -  (3*6),  we  represent  (for  the  continuum) 


and 

En°  =  p2/2m.  (h. 2) 

(These  quantities  are  not  expressed  in  atomic  units.)  In  addition,  we  take 
(19) 


for  the  continuum 
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AXh-H)  1<?2)  '  e1P'r2/K[-9-(X 


with 


»  *  ^  -  p2/2*> 


•  ?• Vv  -  Is  v  27}  #'<  *  r3  -  p2/2"> 

^(p-v)2^  +  ww]  ( \ +  r  -  p2/2m) 


2  (p*7V)2  (Xv  +  “•  -  y/ 2m) 

8ni  n  r3 


(1».3) 


With  this  approximation  the  continuum  contribution  to  Eqs.  (2.1)4)  -  (2.15)  can 
be  evaluated.  For  this  purpose,  we  note  that  the  resulting  approximation  to 
the  first  series  of  both  Eqs.  ( 2 -lij )  and  (2.1f>),  apart  from  a  factor  of  5, 
then  correspond  to  the  expressions  which  are  obtained  for  a  hydrogen  atom 
centered  about  nucleus-3  in  Fig.  1.  These  have  been  given  in  ref.  (9)  and 
will  only  be  summarized  here.  We  obtain  (in  atomic  units) 

,5/2 


M 


(1) 


37C 


1  ( 1  -0-3/2], 

6hrh  r  M  J 


a. 10 


which  can  be  evaluated  in  closed  form,  and 

m  23/2, 

0>c  -  -  ®c(1)^«  -  -1 ST2  \  rl(1-2^»  rl>  '  r0(1-2^  ^  •  tt-5) 

where  rQ  ^  r^  are  both  determined  by  the  condition  that 
r(l  “^Mr)3  * 


(I4 .6) 


The  approximation  for  the  continuum  contribution  to  the  second  series  of 
Eqs.  (2.1)4)  -  (2.15)  entails  the  multiplication  of  Eq.  (it.3)  by  e r3^ 
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by  Eq.  (3*7),  and  an  integration  of  the  product  over  available  momenta.  This 
has  been  carried  out  but,  in  the  interest  of  brevity,  the  details  will  be 
omitted.  We  have  obtained 


M 


c 


(2) 


2  2k+1 
^  k=o  (2k+2)iP(|+k+l) 


,2k+2 


d*M 


2k+2 


(V  k 


), 


(U.7) 


where  R  is  the  internuclear  distance  (in  atomic  units).  With  the  aid  of 
Eq.  (2*7),  we  can  obtain 


(2) 

M 


(M) 


00 

The  continuum  approximations  to  the  normalization  and  energy  are  to  be 
halved  and  added  to  the  discrete  contribution  according  to  Eqs.  (2.II4)  and 
(2.15)*  That  is,  we  take 


(US) 


and 


U.10) 


3>.  Computation 

The  evaluation  of  Eqs.  ( J4 .9)  and  (I4.IO)  were  carried  out  numerically  as 
follows.  For  the  integrals  in  Eqs.  (3.10)  and  (3.1l),  the  integrand  could  be 
expressed  in  terms  of  Bessel  functions,  following  Coulson^^,  and  reduced  to 
radi  al  integrations  on  r, •  These  integrals  were  evaluated  analytically  and 

V.  -> 

(  21) 

expressed  in  terms  of  the  upper  limit  given  by  Eq.  (3.13) v  .  Different  in¬ 
tegrals  arising  from  the  use  of  different  wave  functions  are  related  by  a 
common  value  of  ^  Hence  the  upper  limits  of  their  respective  ranges  of 


13 


integration  are  simply  determined  and  depend,  for  a  given  value  of  ^  only 
upon  the  principal  quantum  number  of  the  wave  function.  The  procedure  followed 
consisted  in  choosing  a  value  of  evaluating  the  appropriate  upper  bounds 

and  substituting  into  the  analytical  expressions  for  the  integrals.  Actually, 
the  discrete  contribution  was  limited  to  states  for  which  n^3.  (Estimates  of 
the  contribution  to  the  normalization  were  made  for  n  =  h  which  led  us  to  be¬ 
lieve  that  their  contribution  was  less  than  0.02.) 

For  a  fixed  value  of  the  quantities  in  Eqs.  (h .It)  -  (lj.8)  were 

evaluated  and  combined  according  to  Eqs.  (U»9)  -  (h.10). 

These  computations  were  carried  out  for  a  fixed  value  of  the  internuclear 
distance,  R,  until  successive  trial  values  of  if  gave  a  value  of  unity  (+0.005) 
for  the  normalization  integral. 

The  procedure  was  repeated  for  different  values  of  the  internuclear  dis¬ 
tance  . 

6.  Results 

The  results  which  have  been  obtained  are  exhibited  compactly  in  Figs. 

(2)  -  (5)«  The  approximations  inherent  in  the  truncation  of  the  discrete 

contribution  and  the  plane-wave  representation  of  the  continuum  contribution 

warrant  no  more  elaborate  presentation.  Nevertheless,  from  Figs.  (2)  and  (3) 

the  computed  values  for  the  potential  energy  of  the  Hg+  molecule,  in  both  the 

ground  gerade  and  unger ade  states,  are  in  good  agreement  with  those  obtained 

( 22) 

from  the  solution  of  the  relevant  Schroedinger's  equation  .  For  the  bound 
state,  a  binding  energy  of  0.085  a.u.  and  an  equilibrium  distance  of  2.10  a.u. 
have  been  obtained;  these  compare  favorably  with  the  quantum  values  of  0.103 
a.u.  and  2.00  a.u.,  respectively. 


l h 


In  Figs,  (h)  and  ($)  are  presented  the  contributions  to  the  normalization 
(of  unity)  made  by  the  various  states.  The  indications  are  that  an  important 
continuum  contribution  is  always  present  for  the  range  of  internuclear  dis¬ 
tances  used  in  these  computations.  Moreover,  its  importance  increases  for 
large  internuclear  distances,  becoming  equal  in  value  to  the  discrete  con¬ 
tribution  when  the  internuclear  distance  becomes  indefinitely  large. 

Although  the  results  which  have  been  obtained  suggest  that  a  sufficiently 
quantum-mechanical  version  of  statistical  theory  can  yield  results  for  mole¬ 
cular  binding  energies  and  bond  distances  which  are  in  good  agreement  with 
those  obtained  from  precise  calculations,  the  method  described  in  this  paper 
holds  little  promise  of  extension  to  polyelectronic  molecules.  In  the  absence 
of  any  self-consistent  field  approximations,  the  resulting  integrals  will  be 
even  more  complicated  than  those  of  the  present  paper.  The  latter  involve 
integrations  over  finite  regions  of  configuration  space  of  one  electron;  in 
polyelectronic  molecules  the  corresponding  regions  for  each  electron  will  be 
correlated  with  those  of  the  others,  rendering  an  analytical  evaluation  of 
the  pertinent  integrals  impractical,  if  not  impossible. 
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